Quantum phase transition in a ID transport model with boson affected hopping: 
Luttinger hquid versus charge-density-wave behavior 
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We solve a very general two-channel fermion-boson model describing charge transport within 
some background medium by means of a refined pseudo-site density matrix renormalization group 
(DMRG) technique. Performing a careful finite-size scaling analysis, we determine the ground-state 
phase diagram and convincingly prove that the model exhibits a metal-insulator quantum phase 
transition for the half-filled band case. In order to characterize the metallic and insulating regimes 
we calculate besides the local particle densities and fermion-boson correlation functions, the kinetic 
energy, the charge structure factor, the Luttinger liquid charge exponent and the single-particle 
excitation gap for a one-dimensional infinite system. 
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The proof of existence of metal-insulator transitions 
(MITs) in generic model Hamiltonians is one of the most 
fundamental problems in solid state theory. While the 
mechanisms that can drive a MIT, such as band struc- 
ture effects disorder Q, Coulomb correlations §], or 
the coupHng to the lattice degrees of freedom [4|, are 
accepted in general, there is only a very small num- 
ber of microscopic models which have rigorously been 
shown to indeed exhibit such a transition. Examples 
are the 3D Anderson tight-binding (disorder) model for 
which an analytical proof of particle localization ex- 
ists [H, or the ID spinless fermion Holstein (electron- 
phonon) model, where the Tomonaga-Luttinger-liquid 
(TLL) charge-density-wave (CDW) MIT has been con- 
firmed numerically by DMRG [6|. Zero-temperature 
MITs triggered by Coulomb interaction are more difficult 
to assess; the prototype half-filled ID Hubbard model, 
e.g., is insulating for all C/ > 0], and only on intro- 
ducing a (particular) long-range hopping the MIT takes 
place at finite interaction strength [8|. In the ID half- 
filled extended Hubbard {U-V) model, there exists at 
most a metallic line at the bond-order- charge-density- 
wave insulator- insulator transition If Coulomb and 
electron-phonon interactions compete, an extended in- 
tervening metallic phase may occur between Mott and 
Peierls insulating states, which allows for a MIT. This 
has been demonstrated for the ID Holstein-Hubbard 
model at half filling [loj . 

Quite recently a novel quantum transport Hamiltonian 
has been proposed , which describes regimes of quasi- 
free, correlation or fluctuation dominated transport. In 
a sense this model parameterizes the correlations inher- 
ent to a fermionic many-particle system, but also the 
couplings to phonon or bath degrees of freedom, by a 
"background medium" that controls the particle's trans- 
port properties. Thus the model captures basic aspects 
of more complicated Hubbard or Holstein Hamiltonians. 
Then it is a legitimate question to ask whether the inter- 
action with the background may even drive a MIT. 



Consider the Hamiltonian [11 1 
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where Hb = ~tbJ2{i,j) fjfii^l + ^j) describes a boson- 
affected nearest-neighbor hopping (oc th) of spinless 
fermionic particles (/j'"^''). A fermion emits (or absorbs) 
a local boson 5] (5,) every time it hops between lattice 
sites This way the particle creates local distortions 
of certain energy in the background. In the case of an 
antiferromagnetic spin background the distortions corre- 
spond to local spin deviations (cf. the motion of a hole 
in the t-J model l^l)- If the background medium is a 



deformable lattice they are basically lattice fluctuations 
(phonons). Other situations such as doped CDWs or ex- 
citon transport in molecular aggregates might be envis- 
aged. In any case the distortions of the background can 
be parameterized as bosons [III ■ The distortions are 
able to relax (compare A with in the t-J model), which 
is described by the second term in ([1]). The third term 
gives the energy of the bosons; the constant energy shift 
NX'^ /ujq guarantees finite energy for — > oo. Performing 
the unitary transformation bi t—^ bi-\- X/uq eliminates the 
boson relaxation term in favor of a second, free-fermion 
hopping channel, H^H = Hi, + Hf+ujo J2i ^l^i ; where 
Hf = —tfJ2{ij) fjfi ^ith tf — 2Xtb/u!o. Hereafter we 
focus on the ID half-filled band case, i.e., fermion number 
Nf/N = 1/2, and take ti, = 1 as energy unit. 

In our model ([T]), the particles have only a charge de- 
gree of freedom. Then, for a tight-binding band struc- 
ture and in the absence of disorder, the formation of a 
CDW is the only possibility for a MIT. The CDW might 
be induced by strong correlations in the background, 
which exist for large luq because (i) coherent transport 
(oc t/) takes place on a strongly reduced energy scale only 
and (ii) incoherent transport is energetically costly 11|. 
Hence, in the limit wq S> 1, an effective Hamiltonian with 
nearest-neighbor fermion repulsion results. By contrast. 
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if the local distortions of the background relax readily 
(i.e., A » 1) and/or the energy of the bosons is small 
(i.e., ujq ^ 1), the free hopping channel can act efficiently 
against any correlation-induced charge ordering. 

Evidence for a MIT comes from a very recent exact 
diagonalization (ED) study of ([T]): Calculating the wave- 
vector-resolved photoemission and inverse photoemission 
spectra the opening of a single-particle excitation gap has 
been observed at Kp — ±tt/2 as A decreases at relatively 



large ujq = 2 [ij]. Of course, dealing with lattices up 
to 16 sites, this does not unambiguously prove the exis- 
tence of a true phase transition which may occur in the 
thermodynamic limit N —^ oo only. 

In this work, we carry out the first large-scale DMRG 
investigation of the two-channel transport model H]). In 
combination with a finite-size scaling analysis this allows 
us to map out the ground-state phase diagram for the 
ID half-filled band case and to characterize the differ- 
ent phases involved. The DMRG is one of the most 
powerful and accurate numerical techniques for studying 
ID fermionic many-body systems [l^. It can be easily 
generalized to treat systems including bosons. Within 
the pseudo-site approach an exact mapping of a boson 
site, containing 2"'' states, to ni, pseudo-sites is per- 
formed [lij . Here we take into account up to rtf, = 4 
pseudo-sites, so that the n^-th local boson density is al- 
ways smaller than 10^^. In addition, we keep m = 1200 
to 2000 density-matrix eigenstates and extrapolate all 
quantities to the to — > oo limit. To test our DMRG im- 
plementation we compared data obtained for small sys- 
tems with previous ED results 1J| and got very good 
agreement: The relative error of the ground-state energy 
\E-EB — E-£)MRg \ was always smaller than 10~^ (for all A at 
ujQ ~ 2); the discarded weight was smaller than 5 x 10^^. 

As indicated by small cluster EDs [l^, at loq — 2.0, 
where fermions and bosons are strongly correlated for 
small A, a MIT might occur in the range of 0.01 < A < 5. 
This is confirmed by DMRG for much larger systems: 
Figs. [1] (a) and (b), showing the variation of the local 
densities of fermions (///i) and bosons {b\h^, respec- 
tively, point towards the existence of a homogeneous 
state (CDW state) for rather large (small) A. Using open 
boundary conditions (OBC), the system is obviously not 
translation invariant, i.e., the local density is inhomoge- 
neous in any case. In the CDW phase, there are two 
degenerate ground states. Within an OBC DMRG cal- 
culation, one of these ground states is picked out by ini- 
tializing the DMRG algorithm, so that the CDW state is 
directly observable in the local density. In the metallic 
regime, on the other hand, the open boundaries reveal 
(strong) Friedel oscillations, which will be algebraically 
reduced, however, as we move towards the interior. Thus 
for large enough system sizes, within the central part of 
the chain, the local density becomes constant (see filled 
symbols in Figs. [T] (a) and (b)) [l7j . 

The CDW structure of the insulating state shows up 
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FIG. 1: (color online) Local densities of fermions {f} fi) (a) 
and bosons {h\hi) (b) for a 128-site system with OBC. Open 
symbols are for A = 0.01 (CDW regime), filled ones for A = 
0.5 (metallic regime). The fermion-boson correlation function 
Xfb{j) is given in panel (c) for a 64-site system with APBC 
(here the discarded weight is 1.4 x 10"^° (7.9 x 10^^°) for 
A = 0.01 (A = 0.50)). In all cases ujq = 2.0. 



also in the fermion-boson correlation function 
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Calculating Xfb{j) at uq = 2.0 for iV = 64 with anti- 
periodic boundary conditions (APBC) [31, we find a 
distinctive alternation for A = 0.01 and a constant value 
away from the "central site" for A = 0.50 which again 
supports the MIT scenario (see Fig. [T] (c)). Note that 
in the latter case there is still a large boson density at 
the particle's nearest-neighbor site, locally enhancing the 
mobility of the carrier. 

Whether the pronounced CDW correlations observed 
for small A and large loq are signatures of true long-range 
order remains an open issue yet. To answer this question, 
we explore the static charge-structure factor. 
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where < q < 27t. If S'c(7r)/iV stays finite in the thermo- 
dynamic limit, CDW long-range order exists. Fig. [2] (a) 
demonstrates that this is the case for A — 0.01, i.e. when 
the distortions of the background relax poorly. By con- 
trast, limAT^oo S{7t)/N = for a = 0.5. This means the 
model ([T]) undergoes a quantum phase transition from 
a metal to an insulator as the relaxation parameter A 
decreases at fixed loq. 
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FIG. 3: (color online) Single-particle gap (squares) and 
TLL parameter Kp as a function of for luq — 2.0 (main 
panel). The inset displays results for smaller ojq and shows 
that (i) no CDW state is found for ujo < uJc and (ii) Kp < 1 
for all Luo, where Kp — > 1 as ljq — > and/or A"^ 0. 



FIG. 2: (color online) Finite-size scaling of several physical 
quantities: (a) static charge structure factor Sc{q) at g = tt, 
(b) kinetic energy parts E^J^, (c) Luttinger liquid parameter 
Kp, and (d) single-particle excitation gap . Data obtained 
for coo = 2 with APBC [(a),(b)] and OBC [(c), (d)] applied. 



Next we investigate the relative importance of the dif- 
ferent transport mechanisms by calculating the kinetic 
energy parts 



(4) 



{\ipo) denotes the ground state). Fig. [5] (b) shows that 

/ 



iJf™ tends to zero at small A, indicating the suppres- 



sion of the coherent transport channel. Note that boson- 
assisted transport is possible for both small and large A 
(^j^km s^a^ys close to —0.6, see panel (b)), and even be- 
comes more pronounced in the CDW phase. 

Finally we determine the TLL charge exponent Kp 
and the single-particle excitation gap (which, for the 
model dJ), equals the charge gap). Kp is proportional 
to the slope of the charge structure factor in the long- 
wavelength limit 5 — > 0+ [igj : 



Ko = TT lim 

9^0 



Sc{q) 



,r 

g = — , N 



(5) 



From this relation we can calculate Kp quite accurately 
using DMRG techniques. As is well-known the ID spin- 
less fermion model with the nearest-neighbor Coulomb 
interaction V at half filling can be mapped onto the ex- 
actly solvable XXZ model. There the TLL charge ex- 
ponent decreases from Kp = 1, as 1^ is enhanced, and 
finally reaches 1/2 at the MIT point We expect 

that this holds also for the ID spinless fermion transport 
model (HI) at half band-filling, even though there are only 
a few analytical or numerical results referring to this for 



coupled fermion-boson systems (for the (Hubbard-) Hol- 
stein model see 0, The single-particle (charge) gap 
can be obtained from 

A,, (N) = EiNf + 1) + E{Nf - 1) - 2E{Nf), (6) 

where E{Nf) and E{NfzLl) are the ground-state energies 
in the Nf- and (TV/ ± l)-particle sectors, respectively, 
with Nf = N/2. Fig. [2] (c) and (d) illustrate the finite- 
size scaling analysis for the TLL parameter (c) and the 
charge gap (d) at wq = 2.0. Both physical quantities 
can be extrapolated by performing a least-squares fit to 
a second-order polynomial in 1/iV. Note that close to 
the MIT points we need larger system sizes, because a 
strong finite-size dependence evolves. In this regime, we 
use chains from = 32 to 256 sites and higher order 
polynomial functions (up to fourth order) to extrapolate 
the data. In doing so, we determine the non-universal 
exponents Kp{X) > 0.5 in the metallic TLL phase where 
Acj = 0, and the finite charge gap A^ > in the CDW 
phase (see panels (c) and (d)). 

In Fig. [3] we display the extrapolated values of the TLL 
exponent and the charge gap as a function of A at fixed 
luq — 2.0. Lowering A, Kp decreases from 1 1/2. 
The point where Kp = 1/2 is reached marks the criti- 
cal couphng for the MIT {X'^itoo = 2) - 20.4). From 
the extrapolated DMRG data it seems that the charge 
gap opens exponentially on entering the insulating phase 
and afterwards rises almost linearly. This result is sim- 
ilar to what is observed for the TLL-CDW transition 
in the anti-adiabatic strong-coupling limit of the spin- 
less fermion Holstein model, which there possesses XXZ- 
model physics (i.e., a Kosterlitz-Thouless transition at 
the spin isotropy point). But note that we find a repul- 
sive particle interaction {Kp < 1) in the metallic phase 
for small boson frequencies as well, i.e. there is no in- 
dication for a pairing instability in the half-filled band 
case. 
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phase typifies a repulsive Luttinger liquid, while the in- 
sulating phase shows CDW long-range order. The phase 
boundary between these states is non-trivial. It would 
be highly desirable to verify the numerical results of this 
paper by an analytical (field theoretical or algebraic) ap- 
proach. 
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FIG. 4: (color online) DMRG phase diagram of the two- 
channel transport model for the ID half-filled band case 
(the dashed line is a guide to the eye). The inset gives the 
phase diagram in the A-uo plane. The MIT point for A = 0, 
tJo(O) ~ 1.38, is obtained from a quadratic fit (dotted line). 



Figure 2] represents the main result of our work: The 
ground-state phase diagram of the ID half-filled fermion- 
boson model (H]) in the A~^-a;,^^ plane. Obviously the 
phase space is divided into two regimes, the metallic TLL 
phase and the insulating CDW phase with long-range or- 
der. We first discuss the limit of large c^o. In this regime 
background fluctuations, which are intimately connected 
with any particle hop, are energetically costly. As a re- 
sult the itinerancy of the particles is suppressed to a large 
extent and charge ordering becomes favorable. Never- 
theless, we find a metallic state even for uq — oo pro- 
vided that A~^ < X~^{u!o ~ oo) (numerically we proved 
the TLL to exist for A^^ < A^^ ~ 6.3 at loq = 1000). 
In this case the system's ability for relaxation (oc A) is 
strong enough to prevent long-range charge order. This 
is reminiscent of the existence of a finite critical cou- 
pling strength gc {g^ = Ep/oJo) in the antiadibatic limit 
(wq oo) of the spinless fermion Holstein model, where 
the TLL phase is realized for g < gd^) iSj- In contrast to 
the TLL-CDW transition in the Holstein model, however, 
the symmetry-broken CDW state is a few-boson state [13] 
(i.e., not a Peierls phase with many phonons involved). 
In the opposite limit of small loo , the background medium 
is easily disturbed by particle motion. Therefore the rate 
of bosonic fluctuations (oc ijJq^) is high. Now we enter 
the fluctuation dominated regime [ll[, and consequently 
CDW order is suppressed. The inset of Fig. |4] shows that 
even for A = 0, i.e. if the explicit A-relaxation channel is 
closed, a metallic state may exist below a finite critical 
energy wo(0). 

To conclude, using an unbiased numerical (DMRG) 
technique, we proved that the very general fermion-boson 
transport model ([T]) displays a correlation-induced metal 
insulator transition at half filling in ID. The metallic 
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